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We measure the polarization ol W bosons from top-quark (t) decays into final states with a 
charged lepton and jets, tt — > W + bW~b — > Ivbqq'b, using the full Run II data set collected by 
the CDF II detector, corresponding to an integrated luminosity of 8.7 fb~ . A model-independent 
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to yield / = 0.726 ± 0.066 (stat) ± 0.067 (syst) and /+ = -0.045 ± 0.044 (stat) ± 0.058 (syst) with 
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assumptions. No significant discrepancies with the standard model are observed. 
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Run I (1992-1996) data sets collected by the CDF and 
DO experiments at Fermilab [l[. Because the mass 
of the top quark is large, beyond-the-standard-modcl 
(BSM) physics contributions can affect the top-quark 
phenomenology in a wide variety of ways: the produc- 
tion mechanisms can be affected, the decay widths can 
be altered, its intrinsic properties can be changed, and 
the experimental signature can be mimicked by a new 
particle of similar mass. Thus a principal goal of the 
Tevatron Run II (2001-2011) program, which produced 
data samples O(100) times larger than Run I, was to 
thoroughly explore the properties of the top quark. This 
Letter reports a measurement of the VF-boson polariza- 
tion from top-quark decay using the full Run II data set 
collected by the CDF II experiment. The focus is on the 
lcpton-plus-jets final state, tt — > W + bW~b — s- Ivbqq'b, 
which provides the most sensitive determination of the 
W polarization due to its high yield, low background, 
and constrained kinematics. At present the most pre- 
cise experimental knowledge of the W^-boson polariza- 
tion comes from the ATLAS experiment [2f and from the 
combination of existing Tevatron results [3(. The results 
reported here supersede the previous CDF measurements 
in Ref. Q and have a total uncertainty comparable to the 
Tevatron combination and a statistical precision a factor 
of 1.6 smaller than Ref. 

The top quark Q almost always decay to a W + boson 
and a b quark [(| via the charged-current weak interac- 
tion whose V — A structure in the standard model (SM) 
specifies the tWb coupling and the resulting VF-boson 
polarization. Due to its large mass, the top quark de- 
cays before hadronizing and thus offers a direct probe 
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of the tWb coupling. At first order in the SM per- 
turbative expansion Q, the W + boson is expected to 
have longitudinal polarization /o = 0.696, left-handed 
polarization /_ = 0.303, and right-handed polarization 
/+ = 3.8 x 10~ 4 for a top-quark mass m t = 172.5 
GeV/c 2 [§], a 6-quark mass m b = 4.79 GcV/c 2 fgj], and a 
IL-boson mass M w = 80.413 GeV/c 2 [§]. Higher-order 
quantum chromodynamic and clcctroweak radiative cor- 
rections, as well as the uncertainties on the values of m t , 
mb, and Mw, change these predictions at the 1-2% rela- 
tive level 0, [l(J ■ The presence of anomalous couplings at 
the tWb vertex, due to contributions from BSM physics, 
can modify the observed W polarization with respect to 
the SM expectations Q. 

In this Letter, three different measurements of the W- 
boson polarization are performed: a modcl-indcpcndcnt 
determination that simultaneously measures fo and /+, a 
measurement of fo for fixed /+=0, which enhances sensi- 
tivity to anomalous tensor couplings, and a measurement 
of f + for fixed fo = 0.70, which enhances sensitivity to 
anomalous right-handed couplings. The analysis assumes 
a top-quark mass of mt = 172.5 GeV/c 2 , consistent with 
the world average value [1] . 

We use a data sample enriched in it — > W + bW~b — > 
Ivbqq'b events, where one of the W bosons decays into 
quark pairs and the other into lepton pairs. The data 
was acquired by the Collider Detector at Fcrmilab (CDF 
II) [ll|], which recorded pp collisions from Fermilab's 
Tevatron operating at y/s = 1.96 TeV. Most of the events 
used in the analysis were collected using inclusive-lepton 
online event selections (triggers) that required a high- 
transverse-momentum (pt) electron or muon in the cen- 
tral (pseudorapidity \rj\ < 1.1) detector region [jjj]. The 
acceptance for it events is increased by also using muon 
events satisfying a trigger that requires large missing 
transverse energy Et [12f with either an energetic elec- 
tromagnetic cluster or two separated jets [l3| ($r+jets 
trigger). After all data quality requirements, the sam- 
ple collected corresponds to an integrated luminosity of 
8.7 fb _1 . 

Candidate events are required to have a single isolated 
electron or muon candidate with Et > 20 GeV; missing 
transverse energy Et > 20 GeV consistent with expecta- 
tions from the undetectable high energy neutrino; and at 
least four energetic jets with Et > 20 GeV and [n| < 2. 
Jets are reconstructed using a cone algorithm 14[ with 
radius Ai? = 0.4 in r) — <fi space, and their energies are 
corrected to the particle level by accounting for detector- 
response non-uniformities as a function of jet r\, for ef- 
fects from multiple pp interactions, and for the hadronic 
jet energy scale of the calorimeter [l5[ . At least one jet 
must be identified as having originated from a b quark 
(b tag) using an algorithm that exploits the long lifetime 
of b hadrons and their large boost from the decay of the 
top quark. We require decay vertices displaced from the 
primary pp- interaction vertex (lil ] . 



The backgrounds to the it signal are from multi-jet 
production (QCD), direct W-boson production in asso- 
ciation with jets (W+jets), and electroweak backgrounds 
(EWK) composed of diboson ( WW, WZ, ZZ) and single 
top-quark production. The IL+jets background includes 
events with 6-tagged 6-quark jets as well as erroneously 
6-taggcd light-flavor or charm-quark jets. The method 
for estimating the background is described in detail in 



Ref. [17| . Table U shows the expected sample compo- 
sition using a it cross section of 7.4 pb from Ref. [18| . 
Events that satisfy one of the high-p^ lepton triggers and 
the ^T+jets trigger are assigned to the lepton-triggered 
sample and removed from the ^?T+jets-triggered sample. 



TABLE I: Number of expected and observed events in 8.7 
fb _1 of data for each of the exclusively-defined samples. A it 
cross section of 7.4 pb is assumed for the signal. 



Process 


Central 


Central 


#r+jets 




e 






it 


923 ± 93 


696 ± 54 


441 ± 44 


VL+jets 


160 ± 41 


125 ± 19 


106 ± 21 


EWK 


36 ± 17 


27 ± 11 


16 ± 8 


QCD 


121 ± 48 


6 ± 2 


7 ± 3 


Total expected 


1239 ± 113 


853 ± 59 


569 ± 50 


Observed 


1226 


804 


544 



The tt signal events are modeled using the her- 
wig [l]| Monte Carlo (MC) generator. The QCD back- 
ground is modeled using data control samples. The ALP- 
GEN 2l|, madevent [22| and MC@NLO 23[ programs, 
with PYTHIA [2(J or herwig supplying the parton-shower 
and fragmentation model, and the full PYTHIA [2(| gen- 



erator, are used to model the remaining backgrounds 
and for estimatin g sy stematic uncertainties. A GEANT- 
based simulation 24j is used to model the response of 
the CDF II detector for these simulated samples. The 
signal and background modeling has been extensively 
checked. The observed data and the predicted signal- 
plus-background distributions for various kinematic vari- 
ables are compared in Fig. [TJ We further validate the 
background model using a high-statistics background- 
dominated data control sample obtained vetoing events 
with a 6-tagged jet. 

To determine the polarization fractions fo, /_, and 
/ + , an unbinned likelihood technique is employed. The 
likelihood is calculated using the theoretical matrix el- 
ements for both the dominant signal process, qq — > it, 
and the main background process, inclusive production 
of IL+jets. The method assumes that pp — > it produc- 
tion is accurately described by the SM and includes the 
physical constraint ' + fi = 1- The technique was 
first developed for measuring the mass of the top quark 
and for determining fo when constraining / + to its SM 
value (2^1 . Wc have extend the technique to enable the 
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• data (8.7 fb ) 
— tf f =0.7, f + =0.0 
■■■ttf =0.7, f t =0.3 
□ Background 
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Di-jet Mass (GeV/c 2 ) 



FIG. 1: The data (points) are compared to the MC predic- 
tion for different W polarization fractions. The background 
predictions are shown as the shaded histograms while the 
signal-plus-background predictions are shown as the open his- 
tograms corresponding to (fo, /+) values of (0.7, 0.0) and 
(0.7, 0.3) for the solid and dashed lines, respectively. The 
four kinematic variables displayed are the leading jet pr, the 
lcpton pr, the invariant mass of the pair of light-quark jets 
from the hadronically decaying W boson, and the cosf9* of 
the leptonically decaying W boson. For the latter two distri- 
butions the jet-parton assignment most consistent with the 
signal hypothesis is shown. 



simultaneous determination of fo and /+ Q ■ The it ma- 
trix element is expressed in terms of the VF-boson polar- 
ization fractions and the cosine of the angle 9* between 
the momentum of the charged lepton or down-type quark 
from the VF-boson decay in the VF-boson rest frame and 
the direction of the top quark. For the signal qq — > it 
process (2(|, the leading-order matrix element is used, 



\M\ 



FeF h (2-p 2 sm 2 6 qt ), 



where g s is the strong coupling constant, 6 q t describes the 
angle between the incoming parton and the top quark in 
the rest frame of the incoming partons, and (3 = v/c 
where v is the velocity of the top quarks in the same 
rest frame. The factors Ft and Fh correspond to top 
quarks with a leptonic and a hadronic VF-boson decay, 
respectively, such that 



Fj 



2irg* v mi 



3m+r, 



^ (2E? + 2>E* h m p 



Id 



(1 



i ) 2 /+ + i (i-cos 2 r)/ 



f(l-/o-/+) 



where gw is the weak coupling constant, mj y is the 
charged lepton-ncutrino invariant mass, Tt is the decay 



width of the top quark, and E£ = - 2r £ — — . The 
hadronic factor Ff, is similar. Since we cannot distin- 
guish between up-type and down-type quark jets from the 
hadronic l^-boson decay, both jet-quark assignments are 
used and an average Fh is calculated from the two per- 
mutations. The background matrix element is approx- 
imated usin g th e sum of VF+jets matrix elements from 
the VECBOS [23 program. 

The polarization fractions are determined by maximiz- 
ing the likelihood function L with respect to fo, /+, and 
the fraction of events consistent with the it signal hy- 
pothesis, C s , 



N 



L(f ,f+,C s \x) = l 



C s 



Ps(xt]fo,f+) 
(A a ( Xi ; fo, f + )) 



+ (1 - C s ) 



Pb{x l ) ' 
(A b {x l ))_ 



where N is the number of observed events, x is the vector 
of observed momenta of the final state partons, and (A s ) 
and (At) are the average acceptances for it and W^+jets 
background events, respectively. The dependence of (A s ) 
on the polarization fractions is properly included. The 
signal probability density, P s , and background probabil- 
ity density, P&, are calculated as described in Ref. [iH 
and integrated over the relevant differential cross section, 
which depends on the matrix elements described above, 
convolved with the proton parton distribution functions 
(PDFs). Poorly-known parton- level quantities are inte- 
grated out. The parton four- momenta arc estimated from 
the measured momenta of the trigger lepton and the four 
highest-transverse-energy Et jets in the event. Detector 
resolution effects are accounted for with transfer func- 
tions derived from simulated it samples. There is an 
ambiguity in the jet-parton assignments and all permu- 
tations are used for each event. When calculating P s we 
fix the top-quark mass to m t = 172.5 GeV/c 2 and scan 
the (fo, f+) parameter space. The calculation of Pb is 
independent of mt, fo, and /+. 

The polarization fractions determined from the likeli- 
hood fit differ from the true polarization fractions be- 
cause the signal and background probability densities 
contain approximations. For example, they do not ac- 
curately account for the effects of extra jets arising from 
initial and final state radiation (ISR/FSR) or for the 
full set of contributing background processes. Thus, the 
results of the likelihood fit arc calibrated with samples 
of it and background events simulated using the sample 
composition of Table [I] and assuming a variety of input 
(fo, f+) values. The mean measured polarization frac- 
tions determined from the simulated samples are plot- 
ted against the true polarization fractions and a cali- 
bration function is determined from a linear fit to the 
resulting curve. For the one-dimensional polarization 
measurements a one-dimensional calibration function is 
employed, while for the simultaneous determination of 
(fo,f+) a two-dimensional calibration function is used. 



G 



The resulting calibration functions are used to estimate 
the true polarization fractions from the measured polar- 
ization fractions determined from the three separate like- 
lihood fits. The uncertainties on the coefficients of the 
calibration functions are included in the method-related 
systematic uncertainties. Even though the likelihood can 
be calculated only for the physical values of /o and /+, 
after calibration the corrected measured values can be 
slightly outside their physical ranges. 

The robustness of the fitting and calibration procedure 
is tested over all physical values of (fo,f+) using sim- 
ulated experiments constructed to have the number of 
observed data events and the sample composition of Ta- 
ble HI No significant biases are observed. However, near 
the physical boundaries the statistical uncertainty is un- 
derestimated by as much as a factor of two. A correction 
to the statistical uncertainty is applied in these regions. 
Assuming the SM, the expected statistical uncertainties 
for the simultaneous measurement, after all corrections, 
are ±0.075 and ±0.047 for f and /+, respectively. 

The sources of systematic uncertainty affecting the 
measurements are summarized in Table [TXJ All system- 
atic uncertainties are determined by performing simu- 
lated experiments in which the systematic parameter in 
question is varied, the default method and calibrations 
are applied, and the shifts in the mean measured po- 
larization fractions are used to quantify the uncertainty. 
All shifts are evaluated using the SM polarization frac- 
tions. The dominant source of systematic uncertainty 
for the simultaneous measurement of /o and / + is due 
to uncertainties on the background shape and normal- 
ization. For the model-dependent measurements several 
sources of systematic uncertainty contribute at a compa- 
rable level. 



TABLE II: Summary of systematic uncertainties. 



Source 


A/o 


A/+ 


A/o 


A/+ 








Simultaneous 


Background model 


0.007 


0.011 


0.049 


0.036 


ISR/FSR 


0.011 


0.017 


0.022 


0.023 


MC generator 


0.012 


0.009 


0.023 


0.011 


Color reconnection 


0.013 


0.010 


0.020 


0.016 


Method- related 


0.014 


0.020 


0.018 


0.016 


Jet energies 


0.016 


0.017 


0.010 


0.022 


PDF 


0.024 


0.013 


0.009 


0.016 


Multiple interactions 


0.009 


0.013 


0.008 


0.014 


Total 


0.040 


0.040 


0.067 


0.058 



The uncertainty on the background model is deter- 
mined by simulating experiments with the mean number 
of total background events increased and decreased by 
one standard deviation while keeping the relative contri- 
butions of the various background sources fixed as given 



in Table |U Half the mean difference between these two 
sets of simulated experiments is assigned as the system- 
atic uncertainty. The background model uncertainty also 
adds in quadrature the largest observed change when 
varying the normalization for each background source in 
turn by one standard deviation, while keeping the to- 
tal background fixed thereby affecting the shape of the 
background distributions. The ISR/FSR uncertainty is 
evaluated using MC samples generated with ISR/FSR 
settings that are amplified or damped relative to the de- 
fault settings. The MC generator uncertainty is evalu- 
ated by comparing between tt MC generated by PYTHIA 
and MC@nlo with parton showering done by herwiG; it 
includes uncertainties from not using the NLO matrix ele- 
ment in the generator, choice of parton shower model and 
modeling of tt spin-correlation. In Ref. [H we only listed 
uncertainty from choice of parton shower model. The 
color reconnection systematic uncertainty [29[ is eval- 
uated using MC samples generated with and without 
color reconnection effects adopting different configura- 
tions [5o| of pythia. The method-related uncertainty 
includes propagating the uncertainty on the fit param- 
eters of the calibration functions, including their corre- 
lations. It also includes the uncertainties related to the 
statistics of the MC samples used to perform the calibra- 
tion. The uncertainties in the jet energy scale corrections 
are propagated through the analysis by varying the cor- 
rections within one standard deviation and recording the 
resulting shifts in the polarization fractions. Variations 
associated with the choice of PDF and their uncertainties 
affect the tt acceptance and are included as a systematic 
uncertainty. The luminosity profile of the MC samples 
does not exactly match that of the data. The associated 
systematic uncertainty is evaluated by varying the MC 
distribution of events containing multiple pp interactions 
so that it matches the data distribution. 

Using the 2574 data events that meet all selection cri- 
teria we perform three measurements of the W boson 
polarization fractions. For the model-independent mea- 
surement we simultaneously determine /o and /+ to be 

fa = 0.726 ± 0.066 (stat) ± 0.067 (syst) 
/+ = -0.045 ± 0.044 (stat) ± 0.058 (syst) 

after all corrections. The correlation between /o and /+ 
is -0.69. 

The two-dimensional likelihood contour obtained from 
the data only includes the statistical uncertainty. The 
final contour, including statistical and systematic uncer- 
tainties, is obtained by analytic convolution of the data 
likelihood with a two-dimensional Gaussian representing 
the systematic uncertainties, resulting in a new likeli- 
hood L syst . Figure [2] shows the point estimate with er- 
ror bars, corresponding to one-dimensional 68.27% confi- 
dence level (C.L.) uncertainties, and the two-dimensional 
68.27% C.L. region, obtained using - ln(L syst /L™f£) = 
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0.5 and 1.15, respectively, where is the maximum 

value for the likelihood L S yst- We estimate a shift of 
T(0.010 ± 0.004) in / and ±(0.012 ± 0.002) in /+ per 
±1 GeV/c 2 shift in the top-quark mass from the central 
value of 172.5 GcV/c 2 . 

For the measurement fixing / + = 0, we obtain after all 
corrections / = 0.683 ± 0.042 (stat) ± 0.040 (syst). For 
the measurement fixing /o = 0.70, we measure after all 
corrections /+ = -0.025 ± 0.024 (stat) ± 0.040 (syst). We 
estimate a shift of ±(0.007 ± 0.002) in / and ±(0.008 ± 
0.001) in f + per ±1 GeV/c 2 shift in the top-quark mass 
from the central value of 172.5 GeV/c 2 . 



0.1 



Data fit: full CDF Run II data set 
■ Best fit and 1 o uncertainty 

Joint 1a coverage in 2-dimensions 
A Standard model Prediction 




FIG. 2: The final results for the simultaneous measure- 
ment of /o and /+ after all corrections and including sta- 
tistical and systematic uncertainties. The square marker 
with error bars shows the point of maximum likelihood 
and the one-dimensional la uncertainties obtained from the 
— ln(Z/ sy st /L^tt ) ~ O-^ con tour. The shaded ellipse corre- 
sponds to the two-dimensional lrr uncertainty obtained from 
the — ln(L syst /L™Jf) = 1.15 contour. The triangle marker 
shows the SM prediction. 

In summary, we present measurements of the polariza- 
tion of W bosons in top-quark decays using the lepton- 
plus-jets final state and the full CDF Run II data sample 
corresponding to an integrated luminosity of 8.7 fb~ . A 
matrix-element technique is used to significantly improve 
the statistical precision relative to previously-used tech- 
niques and is extended to allow for a simultaneous deter- 
mination of /o and /+ in a model-independent manner. 
This result improves the statistical precision on both the 
model-independent and model-dependent determinations 
of /o and /+ by a factor of 1.6 compared to the previ- 
ous CDF measurement Q in the lepton±jets channel. 
Our result is the first VF-polarization measurement us- 
ing the full data set from Tevatron Run II, and is the 
most precise single-channel measurement to date from 
the Tevatron. The results from the model-independent 
and model-dependent measurements of W polarization 
are limited by the size of the systematic uncertainties, 
and have a precision comparable to the combination re- 



ported in Ref. |3j . All the results are consistent with the 
SM. 
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